Proportional Fair Coding for Wireless Mesh 

Networks 

K. Premkumar, Xiaomin Chen, and Douglas J. Leith 



Abstract — We consider multi-hop wireless networks carrying 
unicast flows for multiple users. Each flow has a specified 
delay deadline, and the lossy wireless links are modelled as 
binary symmetric channels (BSCs). Since transmission time, also 
called airtime, on the links is shared amongst flows, increasing 
the airtime for one flow comes at the cost of reducing the 
airtime available to other flows sharing the same link. We 
derive the joint allocation of flow airtimes and coding rates that 
achieves the proportionally fair throughput allocation. This utility 
optimisation problem is non-convex, and one of the technical 
contributions of this paper is to show that the proportional 
fair utility optimisation can nevertheless be decomposed into 
a sequence of convex optimisation problems. The solution to 
this sequence of convex problems is the unique solution to the 
original non-convex optimisation. Surprisingly, this solution can 
be written in an explicit form that yields considerable insight 
into the nature of the proportional fair joint airtime/coding rate 
allocation. To our knowledge, this is the first time that the utility 
fair joint allocation of airtime/coding rate has been analysed, 
and also, one of the first times that utility fairness with delay 
deadlines has been considered. 

Index Terms — Binary symmetric channels, code rate selection, 
cross-layer optimisation, optimal packet size, network utility 
maximisation, resource allocation, scheduling 

I. Introduction 

In this paper, we consider wireless mesh networks with 
lossy links and flow delay deadlines. Packets which are 
decoded after a delay deadline are treated as losses. We derive 
the joint allocation of flow airtimes and coding rates that 
achieves the proportionally fair throughput allocation. To our 
knowledge, this is the first time that the utility fair joint 
allocation of airtime/coding rate has been analysed, and also, 
one of the first times that utility fairness with delay deadlines 
has been considered (also, see [1], [2]). 

In the special cases where all links in a network are loss- 
free or all flow delay deadlines are infinite, we show that 
the proportionally fair utility optimisation decomposes into 
decoupled airtime and coding rate allocation tasks. That is, 
a layered approach that separates MAC scheduling and packet 
coding rate selection is optimal. This corresponds to the 
current practice, and these tasks can be solved separately using 
a wealth of classical techniques. 

However, we show that no such decomposition occurs when 
one or more links are lossy or, one or more flows have finite 
delay deadlines. Instead, in such cases, it is necessary to 
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Fig. 1. An illustration of a single cell wireless LAN with 3 flows. The 
central node represents the Access Point (AP), and the other nodes represent 
the wireless STAtions. 



jointly optimise the flow airtimes and coding rates. Further, we 
show that the resulting allocation of airtime and coding rates 
is qualitatively different from classical results. For example, 
consider a single hop wireless network carrying three flows, 
see Figure 1. Flow fx is a delay-sensitive flow (e.g. video) 
while flows f 2 and / 3 are delay-insensitive flows (e.g. TCP 
data). Transmissions are scheduled in a TDMA manner, and 
the delay deadline for flow fx is one schedule period, while 
the delay deadline for flows f 2 and / 3 is infinite. The channel 
symbol error rate is 10 -2 for all flows, and flows use MDS 
codes for error correction. The proportionally fair airtime 
and coding rate allocation that we show in this paper (see 
Eqns. (21) and (22)) results in the allocation of 41% of the 
airtime for flow fx while flow f 2 and flow / 3 each receive 
29.5%. Observe that the proportionally fair allocation assigns 
unequal airtimes to the flows, which is a notable departure 
from the usual equal-airtime property of the proportional fair 
allocation when selection of delay deadlines and coding rate 
are not included, e.g. see [3]. The optimal coding rate is 0.62 
for flow fx and 0.97 for flows f 2 and f^. The coding rate 
for flow fi is much lower than for flows f 2 and / 3 , since 
a smaller block size must be used by flow fx (and more 
redundant symbols for error-recovery) in order to respect the 
delay deadline. Due to the delay deadline, these optimal coding 
rates yield non-zero loss rates. For flow fx, the packet loss 
rate at the receiver, after decoding, is 20%, whereas flows f 2 
and / 3 are loss free. This highlights an important feature of the 
joint airtime and coding rate utility optimisation. Namely, that 
it allows the throughput/loss/delay trade-off amongst flows 
sharing network resources to be performed in a principled, 
fair manner. Without consideration of coding rate, the trade- 
off between throughput and loss cannot be fully understood 
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or optimally managed. Without consideration of airtime, the 
contention between flows for shared network resources cannot 
be fully captured. 

Proportional fairness can be formulated as a utility maximi- 
sation task, with the utility being the sum of log flow rates. 
Figure 1(b) compares the optimal network utility with that 
obtained with a classical type of approach where all flows are 
allocated equal airtime and the coding rates are chosen based 
on the channel error probabilities alone (this corresponds to 
ignoring the delay deadline of flow /i). It can be seen that 
the optimal approach that we present in this paper potentially 
offers significant performance benefits over classical methods. 

We note that one of the reasons why the joint selection of 
airtime/coding rate has not been previously studied is that the 
proportional fair utility optimisation is non-convex, and hence, 
powerful tools from convex optimisation cannot be applied 
directly. Also, the study of the throughput performance by 
jointly considering the coding and the MAC has not been 
performed before. One of the technical contributions of this 
paper is to show that the proportional fair utility optimisation 
can nevertheless be decomposed into a sequence of convex 
optimisation problems. The solution to this sequence of convex 
problems is the unique solution to the original non-convex 
optimisation. Moreover, this solution can be written in an 
explicit form thereby yielding considerable insight into the 
nature of the proportional fair airtime/coding rate allocation. 

The rest of the paper is organised as follows. The related 
literature on utility optimal resource allocation is discussed in 
Section II. Section III defines the network model; in particular, 
we describe the mesh network architecture, the traffic model, 
and the channel model. We also discuss the transmission 
scheduling model, decoding delay deadline, and the network 
constraints. In Section IV, we obtain a measure for the end- 
to-end packet decoding error, and describe the throughput of 
the network. In Section V, we formulate a network utility max- 
imisation problem subject to constraints on the transmission 
schedule lengths, and discuss the optimization framework. In 
Section VI, we discuss two special cases of networks: delay- 
insensitive and loss-free networks, and show that the tasks of 
obtaining optimal airtimes and coding rates decouple in these 
special cases. We discuss the optimal airtime/coding solution 
with some examples in Section VII. Finally we conclude 
in Section VIII. The proofs of Lemmas and Theorems are 
provided in the Appendix. 

II. Related Work 

We consider a multi-hop Network Utility Maximisation 
(NUM) problem with deadline constraints and with a practical 
model for the PHY layer. By means of channel coding, we try 
to recover a packet from the channel errors. Having a low 
coding rate helps in recovering the packets, but at the cost of 
a small fraction of payload, and at the cost of the transmission 
airtimes of other flows. Thus, we consider the problem of 
resource allocation that answers the following question: how 
to allocate throughput across competing flows with each flow 
seeing different channel conditions and respecting the delay 
deadline. 



The problem of Network Utility Maximisation (NUM) has 
been studied in various contexts, with NUM as a network 
layering tool introduced in [4]. 

Much of the work on NUM is concerned with the flow 
scheduling and throughput allocation that achieves the network 
stability region. This work focuses on throughput and largely 
ignores delay constraints. Resource allocation problems from 
the viewpoint of network control and stability is studied by 
Georgiadis et al. in [5]. Network flow scheduling problems 
are studied in a utility optimal framework by Shakkottai and 
Srikant in [6]. In all these works and the references therein, the 
emphasis is on the MAC layers and above. In [5], an energy 
optimal scheduling problem is studied in which the PHYsical 
layer is also considered. 

Some recent work explicitly includes delay constraints in 
the utility optimisation. In [1], Li and Eryilmaz studied the 
problem of end-to-end delay constrained scheduling in multi- 
hop networks. They propose algorithms based on Lyapunov 
drift minimisation and pricing, and show that by dynamically 
selecting service disciplines, the proposed algorithms sig- 
nificantly outperform existing throughput-optimal scheduling 
algorithms. In [2], Jaramillo and Srikant studied a resource al- 
location problem in ad hoc networks with elastic and inelastic 
traffic with deadlines for packet reception, and obtained joint 
congestion control and scheduling algorithm that maximises 
a network utility. In this work the focus is on congestion 
control and scheduling, with the PHYsical layer considered 
to be error-free. 

A short, preliminary version of the work in the current paper 
was presented in [7]. 

III. Network Model 

A. Cellular Mesh Architecture 

We consider networks consisting of a set of C > 1 cells, 
C = {1, 2, • • • , C} which define the "interference domains" in 
the network. We allow intra-cell interference (i.e transmissions 
by nodes within the same cell interfere) but assume that 
there is no inter-cell interference. This captures, for example, 
common network architectures where nodes within a given cell 
use the same radio channel while neighbouring cells using 
orthogonal radio channels. Within each cell, any two nodes 
are within the decoding range of each other, and hence, can 
communicate with each other. The cells are interconnected us- 
ing multi-radio bridging nodes to create a multi-hop wireless 
network. A multi-radio bridging node i connecting the set of 
cells B(i) = {ci, .., c„} c C can be thought of as a set of n 
single radio nodes, one in each cell, interconnected by a high- 
speed, loss-free wired backplane. See, for example, Figure 2. 



B. Unicast Flows 

Data is transmitted across this multi-hop network as a set T 
= {1, 2, • • • , F}, F > 1 of unicast flows. The route of each 
flow / G J 7 is given by C f = {ci(/) ; c 2 (/),•• • ,c tf (f)}, 
where the source node s(f) e Ci(/) and the destination node 
d{f) G Q / (/)- We assume loop-free flows (i.e., no two cells 
in Cf are same). 
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flow f3 

Fig. 2. An illustration of a wireless mesh network with 4 cells. Cells a, 
b, c, and d use orthogonal channels CHi, CH2, CH3, and CH4 respectively. 
Nodes 3, 5, and 6 are bridge nodes. The bridge node 3 (resp. 5 and 6) is 
provided a time slice of each of the channels CHi & CH2 (resp. CH2 & 
CH 4 for node 5 and CH 2 & CH3& CH 4 for node 6). Three flows /i,/a, 
and f-j are considered. In this example, Cf 1 = {a, b}, Cf 2 = {d, b, a}, and 
C /3 ={ C ,d}. 



constraint that J2f-cec f ^/> c — ^ c wnere T c is the period of the 
schedule in cell c. We consider a periodic scheduling strategy 
in which, in each cell c, service is given to the flows in a 
round robin fashion, and that each flow / in cell c gets a time 
slice of T/.c units in every schedule. 

E. Flow Decoding Delay Deadline 

At the source node s(f) for flow /, we assume that kt 
symbols arrive in each time slot, which allows us to sim- 
plify the analysis by ignoring queueing. Information symbols 
are formed into blocks of Dfkt symbols, where Df € 
{1, 2, 3, • • • } is the number of time slots that the block may 
span. Each block of Dfkf information symbols is encoded 
into a block of DfUf coded symbols, where rif — kf/rf 
symbols, with coding rate < r/ < 1. Here, U/ is the 
number of encoded symbols transmitted in one slot i.e. the 
transmitted packet size. The code employed for encoding is 
discussed in Section IV. The quantity Df is a user or operator 
supplied quality of service parameter. It specifies the decoding 
delay deadline for flow /, since after the flow destination has 
collected at most Df successive coded packets it must attempt 
to decode the encoded information symbols. 



C. Binary Symmetric Channels 

We associate a binary random variable Ef c [b] with the 6'th 
bit transmitted by flow / in cell c. Et c [b] = indicates that 
the bit is received correctly, and Ef tC [k] = 1 indicates that 
the bit is received incorrectly, i.e., the bit is "flipped". We 
assume that J5/ C [l], Ef jC [2], ■ ■ ■ are independent and identi- 
cally distributed' (iid), and ¥{E ftC [b] = 1} = a f . c € [0,0.5). 
That is, we have a binary symmetric channel with cross-over 
probability ct/, c - A transmitted bit may be "flipped" multiple 
times as it travels along the route of flow /, and is received 
incorrectly at the flow destination only if there is an odd 
number of such flips. The end-to-end cross-over probability 
along the route of flow / is therefore given by 



E 



n 



(l-a/.c) 1 -*" 



{x a e{0,l},ceCj: £ x <= is odd} cGC/ 



Note that we can accommodate transmission of symbols from 
any 2 m = M-ary alphabet (i.e. not just transmission of binary 
symbols) by associating m channel uses of the BSC for every 
transmitted symbol. The symbol error probability (for any 
m > 1) is then given by f3f = 1 — (1 — at)" 1 . 

In this channel model, the channel processes across time 
are independent copies of the BSCs. In practice this can be 
realised by means of an interleaver of sufficient depth (after 
the channel encoder), which randomly shuffles the encoded 
symbols, combined with a de-interleaver (before the channel 
decoder) at the receiver. This interleaving and de-interleaving 
randomly mixes any channel fades, which can then be mod- 
elled as independent channel processes across time. 

D. Flow Transmission Scheduling 

A scheduler assigns a time slice of duration Tf >c > time 
units to each flow / that flows through cell c, subject to the 



F. Network Constraints on Coding Rate 

For flow / in cell c, let Wf, c be the rate of transmission 
in symbols/second, which is determined by the modulation 
and spectral bandwidth used for signal transmission and the 
within-cell FEC used. Each cell c € Cf along the route of flow 
f allocates an airtime of at least in order to transmit the 

w f,c 

packets of flow /. Let T c := {/ G J- : c E Cf} be the 
set of flows that are routed through cell c. We recall that the 
transmissions in any cell c are scheduled in a TDMA fashion, 
and hence, the total time required for transmitting packets for 
all flows in cell c is given by J2fer Since, for cell c, the 
transmission schedule interval is T c units of time, the encoded 
packet size nf must satisfy the schedulability constraint 

/Ik w f- 



Note that since we require sufficient transmit time at each cell 
along route Cf to allow rif coded symbols to be transmitted in 
every schedule period, there is no queueing at the cells along 
the route of a flow. 

IV. Packet Error Probability 

Each transmitted symbol of flow / reaches the destination 
node erroneously with probability f3f. Hence, to help protect 
against errors when recovering the information symbols, we 
encode information symbols at the source nodes using a block 
code (we note here that a convolutional code with zero- 
padding is also a block code). An (n,k,d) block code has 
the following properties. The encoder takes a sequence of 
k information symbols as input, and generates a sequence 
of n > ft coded symbols as output. The decoder takes a 
sequence of n coded symbols as input, and outputs a sequence 
of k information symbols. These information symbols will 
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be error-free provided no more than L^y^J °f tne coded 
symbols are corrupted. The Singleton bound [8] tells us that 
d ^ n—k + 1, with equality for maximum-distance separable 
(MDS) codes. Thus, an MDS code can correct up to 



(1) 



d-1 




n — k 


2 
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errors. Examples for MDS codes include Reed-Solomon codes 
[8], and MDS-convolutional codes [9]. In [9], the authors 
show the existence of MDS-convolutional codes for any 
code rate. Hereafter, we will make use of Eqn. (1), and so, 
confine consideration to MDS codes. However, the analysis 
can be readily extended to other types of code provided a 
corresponding bound on d is available. 

Consider a coded block of flow / and let i e 
{1, 2, • • • ,Dfiif} index the symbols in the block. Let Ef[i] be 
a binary random variable which equals when the i'th coded 
symbol is received correctly and which equals 1 when it is 
received corrupted. F{E f [i] = 1} = /3 f and F{E f [i] = 0} = 
1 — f3f. From Eqn. (1), the probability of the block being 
decoded incorrectly is given by 



Dfn 



E E 



i > 



D f n f - D i k f 



The symbol errors Ef[l], Ef[2], ■■■ , E^Djrif] are i.i.d. 
Bernoulli random variables, and so, the J2iJi f -^/W * s a bi~ 
nomial random variable. Hence, the probability of a decoding 
error can be computed exactly. However, the exact expression 
is combinatorial in nature, and is not tractable for further 
analysis. We therefore proceed by obtaining upper and lower 
bounds on the error probability, and show that the bounds are 
the same up to a prefactor, and that the prefactor decreases 
as the block size DfUf increases. Hence, we pose the NUM 
based on the upper bound on the error probability. Also, we 
relax the following constraints: rif e Z + and kf € Z +) and 
allow them to take positive real values, i.e., n/ e K + and 
k f e M+. 

A. Upper and Lower Bounds 

Lemma 1 (Upper Bound). The end-to-end probabUity e/ of 
a decoding error for flow f satisfies 

e f < exp(-Dfn f I Ef[1] (xf,6f)) (2) 
=: e f (6f,n f ,r f ). 

where Xf := ^—^ L , ff = kf/rif is the coding rate, Of > is 
the Chernoff-bound parameter and the function lz(x;0) := 
Ox — ln(E [e ez ] ) is called the rate function in large deviations 
theory. 

Proof: See Appendix A. ■ 

Lemma 2 (Lower Bound). The end-to-end probability e/ of 
a decoding error for flow f satisfies 



where 



e f >Te^(-DfnfI(B(x f )\\B(Pf))) 



(3) 



r = 



exp(-D f n f H(B(xf))) 



and Xf :— 



l—re 



B(x) is the Bernoulli distribution with 
parameter x, H(V) is the entropy of probability mass function 
(pmf) V, and I(V\\Q) is the information divergence between 
the pmfs V and Q. 

Proof: See Appendix B. ■ 

B. Tightness of Bounds 
It can be verified that 

I Ef[ i] (xf,9f)=efXf-]n(l-Pf+P f e e f) 

Since Of > is a free parameter, we can select the value that 
maximises Ie s [i] { x f\0f) and so provides the tightest upper 
bound. It can be verified (e.g. by inspection of the second 
derivative) that Ie,[i] ( x f>@f) is concave in Of and so the 
KKT conditions are necessary and sufficient for an optimum. 
The KKT condition here is 

dI Ef[1] (x f ;0f) _ f3 f e e f 



86 f 

which is solved by 



l-Pf + Pfe* 



= 



7 f 



= In 



In 



1 



x f 



provided Xf > (3f. Substituting for 0j, 

max/ B/ [i] {xf,0f) = I Ef[1] (xf,e* f ) 



Xfln(*l) + (1- X f)l*(\ 
I{B{xf)\\B(i3f)) 



1-Xf 



and by Lemmas 1 and 2, the probability e/ of a decoding 
error satisfies 



Te -D,n s I(B{xf)\\B(f) s )) < ~ < e - 



D f n f I(B(x f )\\B(i3,)) 



It can be seen that the upper and lower bounds are the same 
to within prefactor T, and the gap between these bounds 
decreases exponentially as the block size DfUf increases. 

V. Network Utility Optimisation 

We are interested in the fair allocation of flow airtimes 
and coding rates amongst flows in the network. Other things 
being equal, we expect that decreasing the coding rate r/ (i.e., 
increasing the number DfUf — Dfkf of redundant symbols 
transmitted) for flow / will decrease the error probability ef, 
and so increases the flow throughput. However, decreasing 
the coding rate increases the coded packet size Dfnf, and 
so increases the airtime used by flow /. Since the network 
capacity is limited and shared by other flows, this generally 
decreases the airtime available to other flows and so decreases 
their throughout. Similarly, increasing the packet size Dfkf 
of flow / increases its throughput but at the cost of increased 
airtime and a reduction in the throughput of other flows. We 
formulate this tradeoff as a utility fair optimisation problem. 
In particular, we focus on the proportional fair allocation since 
it is of wide interest and, as we will see, is tractable, despite 
the non-convex nature of the optimisation. 
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The utility fair optimisation problem is 



max 
0,(n,x) 



U(d,(n,x)) 



subject to < T r . 



(4) 

Vc e C (5) 



f:cec f 

Of > 0, V/ e T (6) 

x f <Xf V/eJ (7) 

Xf > Xf V/e7 (8) 

with 6 := [#/]/e;r the vector of Chernoff parameters, n := 
[nf]f<zjr the vector of flow packet sizes, and x := [xf}f<zjr 
the vector of flow coding rates (where we recall that Xf = 
(1— r/)/2). Eqn. (5) enforces the network capacity (or the flow 
schedulability) constraints, Eqn. (6) the positivity constraint on 
the Chernoff parameters, and the constraints Eqns. (7)-(8) are 
introduced for technical reasons that will be discussed in more 
detail shortly. 

For proportional fairness, we select the sum of the log of 
the flow throughputs as our network utility U. For flow / 
the expected throughput is Dfkf(l — e/) symbols in every 
time interval of duration DfT d ^ (we recall that d(f ) is the 
destination cell of flow /), which is the same as kf(l — ef) 
symbols every time interval of duration T d /f\, where Dfkf is 
the information packet size and ef the packet decoding error 
probability. As the exact expression of e/ is intractable, we 
use the upper bound for ef which is e/. Thus, the objective 
function is given by 

U (6, (n, x)) := ]T In (k f (1 - e f (6 f ,n f , x f ))) 

= J2 ln ( n f r f _ e f( e f> n f> x f))) 

= ^ ln(n/(l - 2a;/) (1 - e f (0 f ,n f ,x f ))) 

= E ln («/) + E ln d-2x / ) 

+ ^ j \n(l-e s {8 s ,n s ,x } )). 

The optimisation problem yields the proportional fair flow 
coding rates and coded packet size n/. Since the PHY trans- 
mission rates Wf iC are known parameters, the coded packet 
size is proportional to the airtime used by a flow (i.e., the 
airtime is given by n//w/ iC ). 

A. Non-Convexity 

The objective function U(9,(n,x)) is separable in 
(6f,(nf,Xf)) for each flow /. However, it can be readily 
verified that ln(l — ef(6f,rif,Xf)) is not jointly concave in 
(Of, (nf,Xf)), and so, the optimisation is non-convex. Hence, 
the network utility maximisation problem defined in Eqns. (4)- 
(8) is not in the standard convex optimisation framework. 

B. Reformulation as Sequential Optimisations 

We proceed by making the following key observation. 



Lemma 3. . For convex sets y and Z, and for a function 
f : y x Z — >• R that is concave in y e y and in z G Z, 
but not jointly in (y, z), the solution to the joint optimisation 
problem 



is unique, 



max f{y,z) 
and is the same as the solution to 



max max f(y, z), 
zez yey 



(9) 



(10) 



iff(y*(z), z) is a concave function of z, where for each z G Z, 

y*(z) := argmax/(y,2:). 
yey 

Proof: See Appendix C. ■ 
This lemma establishes conditions under which we can 
transform a non-convex optimisation into a sequence of con- 
vex optimisations. Roughly speaking, we proceed by optimis- 
ing over each variable in turn and substituting the optimal 
variable value that is found back into the objective function. 
This creates a sequence of objective functions. Provided each 
member of this sequence is concave in the variable being 
optimised (but not necessary jointly concave in all variables), 
the solution to the sequence of convex optimisations coincides 
with the solution to the original non-convex optimisation. Ev- 
idently, the condition that concavity holds for every objective 
function in this sequence is extremely strong. Remarkably, 
however, we show that it is satisfied in our present network 
utility optimisation. 

C. Optimal 6* (x f) 

Taking a sequential optimisation approach, we begin by first 
solving the optimisation 

max U (6, (n, a;)) 

subject to 6 f > 0, V/ € T 

given packet sizes n e and coding rates x e [A/,A/] F . 
The objective function is separable and concave in the OfS. 
The partial derivative of U (0, (n, x)) with respect to Of is 
given by 



dU{6,{n,x)) efDfnf 



80 f 



1 - 



e f 



x f - 



l-0 f +p f e°t 



(11) 



Setting this derivative equal to zero, provided Xf > f3f this is 
solved by 



(3 f e 



or, e e *f 



or,0*f 



x f 

Xf 1 - Pf 
Pfl-Xf 



In 



I-/?/ 



(12) 



Observe that in fact 0*j is function only of Xf and not both 
nf and Xf. The requirement for Xf > f3f ensures that 0j- > 0. 
When Xf < (if, the derivative (Eqn. (11)) is negative for all 
Of > 0. In this case, the optimum Of is zero which yields 
an error probability e/ of one. Thus, for error recovery we 
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require Xf > /3/ i.e. the coding rate 77 < 1 — 2/3/, and for a 
non-empty feasible region in the NUM problem formulation 
in Eqns. (4)-(8) the constraints on Xf should satisfy the 
following: A/ < 1 — 2/3/ and A/ > 0. We note that the capacity 
region for a BSC having a cross-over probability a / with an 
m-ary signalling is (0, m(l — H(a.f))), and the coding rate 
1 — 2/3/ lies in the capacity region. 



D. Optimal (n f ,x f ) 

The next step in our sequential optimisation approach is to 
solve 

max U (0*(x), (n, x)) 

71/ 



subject to 



E 

f-cec f 

x f < A/ 
x f > X f 



< T„. 



W f,c 



Vce C 

V/eJ 
V/eJ 



That is, we substitute into the objective function for the 
optimal 6*f found in Section V-C. Defining // = Xf In (jjj^) + 

(l-.x/)ln(^), 

U (6* (x), (n, x)) = 2 In (n/) + £ In (1 - 2x f ) 



fer 



lnfl 



-DfUflf^ 



It can be verified that [/ (0*(a;), (n, a?)) is not jointly concave 
in (n,x). To proceed, we therefore rewrite the objective in 
terms of the log-transformed variables n/ = ln(n/) and 
If = ln(7/). Observe that the mapping from rtj to n/ is 
invertible and similarly the mapping from Xf to //. Since // 
is a monotone increasing function of Xf (this can be verified 
by inspection of the first derivative), the inverse mapping 
from If to Xf exists and is one-to-one. With the obvious 
abuse of notation, we write inverse map as Xf(If). In terms 
of these log-transformed co-ordinates, the objective function 
is U (6*(I), (tz, ifj. We note that the problem defined in 
Eqns. (4)-(8) is equivalent to the problem, 



u(e,(n,i)) 

subject to ^ 



max 

0,(n,I) 



f-cec f 



W f,c 



Of > 0, 

if < h 
If > x f 



VceC 



(13) 

(14) 

(15) 
(16) 



and hence, by Lemma 3, the solution to the log-transformed 
problem is the same as that of the problem defined in 
Eqns. (4)-(8). We solve the maximisation problem by convex 
optimisation method. We show that the objective function is 
jointly concave in (n,l) in the following Lemma. 



Lemma 4. 

U (e*(I), (n, /)) = £ hf + £ In (l - 2x f (If)) 



+ 



is jointly concave in rif and If. 



Proof: See Appendix D. ■ 
Hence, we have the following convex optimisation problem 

max U(e*(l),(nj)) (17) 

(nj) V ' 



Eg™/ 
< T c 

f.cec f Wf ' c " 

If < A/ 
If>X f 



VceC 



(18) 



V/ e T (19) 
V/ G T (20) 



We solve the above maximisation problem using the La- 
grangian relaxation approach. The Lagrangian function of the 
problem is given by 



L (n,I,p,u, uj 



nf 



^EW/).^y)-& E Wf 

feJ 7 cec \/e^ c 

fer fer 

where p > 0, v > 0, and v > are Lagrangian multipliers 
corresponding to the constraints given in Eqns. (18)-(20). The 
channel error probabilities /3/s are strictly positive, and the 
channel coding rates are always assumed to be in the interior 
of the feasibility region. Hence, the constraints for the channel 
coding rate given in Eqns. (19), (20) are not active at the 
optimal point, and the Lagrangian costs VfS and VfS are zero. 
Thus, the shadow costs corresponding to these constraints will 
not appear in the Lagrangian relaxation. 

Since the optimisation problem falls within convex optimi- 
sation framework, and the Slater condition is satisfied, strong 
duality holds. Hence, the KKT conditions are necessary and 
sufficient for optimality. Differentiating the Lagrangian with 
respect to n/ at n/ = n*j, and setting equal to zero yields the 
KKT condition 



E 

cec f 



Dfnllfjxye-W'W p c n* f 



fn*I f (xj) 



ceCf 



w f,. 



(21) 



Similarly, the KKT condition for Ij is 

D/exp(n} + /*) e -D/exp(n}+7;) 



1 - 2x* f 9* f 



or, 



I _ e - D f cxp(n*+/p 

D fn y f ix})e- D ' n VW 2 If{x* s ) 



? f n* f I f (x* f ) 



1 - 2x* f 0*f(x* f ) 



(22) 
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Combining eqns. (21) and (22), yields 



and the projected subgradient descent update is 



E 

cec f 



Peri} 

W f,c 



1 - 2x* 0*{x*)' 



(23) 



Observe that the LHS is a function of and the RHS is a 
function of x*f. Thus, the choice of packet size parameter 
and coding rate parameter x*j are in general coupled. 



E. Distributed Algorithm for Solving Optimisation 

Given the values of the Lagrange multipliers p* , the solution 
to Eqn. (23) specifies the optimal packet size and coding rate. 
To complete the solution to the optimisation it therefore re- 
mains to calculate the multipliers p* . These cannot be obtained 
in closed form since their values reflect the network topology 
and details of flow routing. However, they can be readily found 
in a distributed manner using a standard subgradient approach. 

We proceed as follows. The dual problem for the primal 
problem defined in Eqn. (17) is given by 



mm 

p>0 



D(P), 



where the dual function D(p) is given by 

D(p) = max £ U f (6* f (I f ),(n f ,I f )) 



(n,I) 



fer< 



W f,c 



(24) 



£ U f {e* f (I f (p)),(n* f (p),I* f (p))) 
e «J(p) 




cec \ feFc 
From Eqn. (24), for any (n, I), 



W f,c 



D(p) > J2Uf(0}(If), (n f ,I f )) +J>c 



cec 



and in particular, the dual function D(p) is greater than that 
for If = If(p) for some arbitrary p, i.e., 

D(p) > J2 Uf (e}(I}(p)),(n}(p)J}(p)) 



5> h-E 



e'V 



= I>(P) + S (Pc " Pc) T c - £ 



;™/(p) 



(25) 



Thus, a sub-gradient of £)(•) at any p is given by the vector 

n* f (p) 



E 



cGC 



Pc(» + 1) = 



PS) - 7 • ( To ~ 



n}(p(i)) 

W f-c 



where 7 > is a sufficiently small stepsize, and [/(a;)] + := 
max{/(cc), 0} ensures that the Lagrange multiplier never goes 
negative (see [10]). 

The subgradient updates can carried out locally by each 
cell c since the update of p c only requires knowledge of the 
packet sizes n* f -(p(i)) of flows / e T c traversing cell c. Thus, 
at the beginning of each iteration i, the flow source nodes 
choose their packet sizes as Dfri*f(p(i)) and the coding rates 
as 1 — 2x*f(p(i)), and each cell computes its cost based on 
the packet sizes (or equivalently the rates) of flows through 
it. The updated costs along the route of each flow are then 
fed back to the source nodes to compute the packet size and 
coding rate for the next iteration. 

Observe that the Lagrange multiplier p c can be interpreted 
as the cost of transmitting traffic through cell c. The amount 
of service time that is available is given by A = T c — 



E 



-. When A is positive and large, then the La- 



grangian cost p c decreases rapidly (because the dual function 
D(-) is convex), and when A is negative, then the Lagrangian 
cost p c increases rapidly to make A > 0. We note that the 
increase or decrease of p c between successive iterations is 
proportional to A, the amount of service time available. Thus, 
the sub-gradient procedure provides a dynamic control scheme 
to balance the network load. 

The resulting distributed implementation of the joint air- 
time/coding rate optimisation task is summarised in Algorithm 
1. 

Algorithm 1 Distributed implementation of joint air- 
time/coding rate optimisation. 



Each cell c runs: 
loop 

1. Pc(* + 1) = 



Pc{i) - 7 ' \ T C 



E 



w f,c 



end loop 

The source for each flow f runs: 
loop 

1. Measure EceC/ w^ 7 tne a S§ re 8 ate cost °f usm g me 
cells along the route of flow /. E.g. if each cell updates 
the header of transmitted packets to reflect this sum, it can 
then be echoed back to the source by the flow destination. 

2. Find the unique packet size n*j and coding rate x*j that 
solve (23). Since there are only two variables, a simple 
numerical search can be used. 

end loop 



VI. Two Special Cases 

A. Delay-Insensitive Networks 

Suppose the delay deadline Df — > 00 for all flows. For 
any positive bounded n^Ip i.e., < n*j.Ij- < 00, the LHS of 
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Eqn. (22) can be written as 

D f n* f r f e- D f n 'f r f 



e 1 i t — 1 

oo (D f n* f Ipi 



oo (D } n*I*.y- 



lim 



-D f n}I} 



1 - e" 



= 0. 



(26) 



Thus, the asymptotic optimal coding rate as the delay 
deadline requirement Df — > oo is the solution to 

_2 Mgf) _ 

i - 2x* e*(x* f ) 

Since /?/ < < 1/2, it is sufficient to find the solution to 



(27) 



Note that 



lim 



0*(a; 



and hence, lim 



V 

* 

x f 



0. 



(28) 



Since this is the limiting solution and x*j > f3f, one can use 
x*j = f3f + e for some arbitrarily small e > 0. Similarly, from 
Eqn. (21), the asymptotic optimal packet size nj- as Df — > oo 
is 

where the multipliers p c are obtained, as before, by subgradient 
descent 

+ (30) 

fer c 



Pc(» + 1) 



W f-c 



Observe that the optimal coding rate x*j — fif + e which is 
given by the solution of Eqn. (28) is determined solely by 
the channel error rate (3f of flow /. It is therefore completely 
independent of the other network properties. In particular, it 
is independent of the packet size n*j used, of the other flows 
sharing the network and of the network topology. Conversely, 
observe that the optimal packet size n*f in Eqn. (29) and 
Eqn. (30) is dependent on the network topology and flow 
routes, but is completely independent of the error rate /?/ 
and coding rate Xf. That is, in delay-insensitive networks, 
the joint airtime/coding rate optimisation task breaks into 
separate optimal airtime allocation and optimal coding rate 
allocation tasks which are completely decoupled. Our optimi- 
sation therefore yields a MAC/PHY layering, whereby airtime 
allocation/transmission scheduling is handled by the MAC 
whereas coding rate selection is handled by the PHY, with no 
cross-layer communication. It is important to note, however, 
that this layering does not occur in networks where one or 
more flows have finite delay-deadlines; see Section VII for a 
more detailed discussion. 



B. Loss-Free Networks 

Suppose the channel symbol error rate /3f = for all flows. 
From Eqn. (12), we observe that 



lim 6* f = oo, 



(31) 



and this yields e/ = for all flows. The objective function 
in Eqn. (17) degenerates to 2~^/e;r m ( n /(l — 2a;/)). We note 
that for any x* f > /?/ > 0, as /?/ I 0, If(x* f ) -> oo. Hence, 
the LHS of Eqn. (22) becomes, 



D f n* f If{x))e- D ' n P'W 

lim J -= — - = , = 0. 



/3j40 



I _ p-Df"-} 1 /^}) 



(32) 



In the same was as in Eqn. (28), this limit can be achieved 
by x*f = (i.e. = 1). Similarly, the optimal packet size is 
-. This optimal packet size is identical to 



that for delay-insensitive networks, see Eqn. (29), and it can be 
verified that in fact it corresponds to the classical proportional 
fair rate allocation for loss-free networks, as expected. 

VII. Examples 

A. Single cell 

We begin by considering network examples consisting of a 
single cell carrying multiple flows. The network topology is 
illustrated schematically in Figure 1 and might correspond, for 
example, to a WLAN. 

1) Mix of delay-sensitive and delay-insensitive flows: 
Suppose the flows in the network belong to two classes, one of 
which is delay-sensitive and has a delay-deadline D whereas 
the other is delay-insensitive i.e. has an infinite delay deadline. 
These classes might correspond, for example, to video and 
data traffic. Figure 3(a) plots the optimal airtime allocation 
as the delay deadline D is varied. In this example, there is 
a single delay-sensitive flow and two delay-insensitive flows, 
and the airtime allocation is shown for the delay-sensitive 
flow and for one of the delay-insensitive flows (both receive 
the same airtime allocation). As expected, it can be seen 
that the airtime allocations of the delay-sensitive and delay- 
insensitive flows approach each other as the delay deadline 
D is increased. However, it is notable that they approach 
each other fairly slowly, and when the delay deadline is low 
the airtime allocated to the delay-sensitive flow is almost 
50% greater than that allocated to a delay-insensitive flow. 
This behaviour is qualitatively different from the classical 
proportional fair allocation neglecting coding rate and delay- 
deadlines, which would allocate equal airtimes to all flows. 
By taking coding rate and delay deadlines into account, our 
approach allows the resource allocation to flows with different 
quality of service requirements to be carried out in a principled 
and fair manner. 

Figure 3(b) plots the optimal airtime allocation as the 
number of delay-insensitive flows is varied. It can be seen 
that the airtime allocated to each flow decreases as N is 
increased, as expected since the number of flows sharing the 
network is increasing. Interestingly, observe that the airtime 
allocated to the delay-sensitive flow is a roughly constant 
margin above that allocated to the delay-insensitive flows. The 
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(a) Optimal airtime allocation vs delay deadline D, N=2. 



(a) Optimal airtime allocation vs channel symbol error rate for 
flow 1, symbol error rate for flow 2 is held fixed at 10 — 2 . 



- Delay-sensitive flow 

- Delay-insensitive flows 




Number N of delay-insensitive flows 

(b) Optimal airtime allocation vs N. D = 1. 

Fig. 3. Single WLAN with one delay-sensitive flow and N delay-insensitive 
flows. Delay sensitive flow has delay deadline D, delay-insensitive flows have 
infinite delay deadlines. Raw channel symbol error rate is 10 — 2 for all flows, 
PHY rate for all flows is 10 symbols per schedule period. Optimal airtimes 
are given as a proportion of the schedule period. 




PHY rate tor flow 1 (symbols/schedule) 



(b) Optimal airtime allocation vs PHY rate of flow 1, PHY rate 
for flow 2 is held fixed at 10 symbols/schedule. 

Fig. 4. Single WLAN with two delay-sensitive flows, both with delay 
deadline D = I. In upper figure, PHY rate for both flows is 10 symbols 
per schedule period and channel symbol error rate for flow 1 is varied. In 
lower figure, channel symbol error rate for both flows is 10 — 2 , and PHY rate 
for flow 1 is varied. 



delay-sensitive flow is therefore "protected" from the delay- 
insensitive flows. However, in contrast to ad hoc approaches, 
this protection is carried out in a principled and fair manner. 

2) Mix of near and far stations: Consider now a situation 
where all flows have the same delay deadline D, but where for 
some flows the sources are located close to the destination and 
for other flows the sources are further away. We therefore have 
two classes of flows, one with a higher channel symbol error 
rate than the other when both use the same PHY rate. Figure 
4(a) plots the optimal airtime allocation for a flow in each 
class as the channel error rate for one class is varied. When the 
channel error rates for both classes is the same (f3f = 10~ 2 ), 
it can be seen that the airtime allocation is the same. As the 
channel error rate decreases, the airtime allocated to flow 1 
decreases. Conversely, as the channel error rate increases, the 
airtime allocated to flow 1 increases. 

Figure 4(b) plots the optimal airtime allocation when flows 
in both classes have the same channel error rate but different 
PHY rates i.e. where the PHY modulation has been adjusted 
to equalise the channel error rates. When the PHY rates 
are the same (w/. c = 10 symbols per schedule period), the 
airtime allocation is the same to both classes. As the PHY 
rate is increased, the airtime allocation for flow 1 decreases. 



Conversely, as the PHY rate is decreased, the airtime allocation 
for flow 1 increases. Again, note that this is qualitatively dif- 
ferent from the classical proportional fair allocation neglecting 
coding rate and delay-deadlines which would allocate equal 
airtimes to all flows. 

3) Unequal Airtimes: The basic observations in these ex- 
amples apply more generally. In particular, as noted above, 
in a loss-free, delay-insensitive single-cell network the pro- 
portional fair allocation is to assign equal air-time to all 
flows ([3] and Section VI-B). However, when delay deadlines 
are introduced and/or links are lossy, we see an interesting 
phenomenon. 

Lemma 5. The optimum rate allocation x* (or equivalently 
r*) is not equivalent to an equal air-time allocation. 

Proof: See Appendix E ■ 
In particular, flows that see a better channel get less air- 
times than flows that see a worse channel. 

B. Multiple cells 

We now consider a mesh network consisting of TV cells 
carrying N + l flows in the well-studied Parking Lot topology. 
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/ flow 1 
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\ flow 2 




flow 3 \ / flow 4 


1 




2 3 



4N+1 



N-l 



Fig. 5. A linear Parking Lot network with TV cells and TV + 1 flows (one 
multi-hop flow and TV single-hop flows). 



The network topology is illustrated in Fig. 5. The flows in this 
network can be assigned to two classes: class 1 consists of the 
TV-hop flow, and class 2 consists of the single-hop flows 2, 
3, ■ • • , N + 1. Each cell has the same schedule period, i.e. 
T c = T,VceC. 

1 ) Impact of number of hops: Suppose that both classes 
of flow use the same symbol transmission PHY rate and 
experience the same loss rate in each cell. Then the TV-hop 
flow will experience a higher end-to-end symbol error rate 
that the single hop flows, and the loss rate will increase with 
N. Fig. 6 plots the ratio of optimal airtimes allocated to each 
class of flow versus TV. Results are shown for three delay 
deadline requirements: both classes of flow are delay-sensitive 
with delay deadline D\ = D% = 1; class 1 is delay-sensitive 
{Di = 1) while class 2 is delay-insensitive (D 2 = 10 5 ); class 

1 is delay-insensitive while class 2 is delay-sensitive. It can 
be seen that in the first case, where both classes have the 
same delay deadline, the ratio of airtimes is larger than 1. 
This is in accordance with the previous observation that flows 
with poorer channel conditions are allocated more airtime than 
flows with better channel conditions. In the second case, where 
class 2 is delay-insensitive (D2 — 10 5 ), additional airtime 
is allocated to class 1, the delay-sensitive flow, which also 
corresponds with the single cell analysis. In the third case, 
where class 1 is delay-insensitive (Di = 10 5 ) and class 2 is 
delay-sensitive, it can be seen that class 2 flows are allocated 
slightly more airtime that the class 1 flow. Interestingly, 
however, observe that the airtime allocated to the class 1 flow 
is insensitive to the number TV of hops. This contrasts with 
the behaviour when the class 1 flow is delay-sensitive. 

2) Impact of different flow PHY rates: Now consider a 
situation where the number of cells TV = 3 and all flows 
have the same delay deadline D\ = D2 = 1. Flow 2 and flow 
4 have symbol error rate 1CP 4 , and flow 1 and flow 3 have 
symbol error rate 2.5 x 10 -1 . We classify the flows into three 
sets: class 1 consists of multi-hop flow 1, class 2 consists of 
single-hop flows 2 and 4, class 3 consists of single-hop flow 
3. Let w denote the PHY rate used used by class 1 and class 

2 flows, and 1/73 denote the PHY rate used by the class 3 flow. 
Fig. 7 plots the optimal coded packet size versus the ratio 
W3/11). We begin by observing that when w 3 /w — 1, all flows 
have the same PHY rate and it can be seen that flows in classes 
2 and 3 are allocated the same packet sizes (and so the same 
airtime). Hence, although the flow in class 3 crosses a much 
more lossy link than the flows in class 2, the optimal allocation 
ensures that all of the single-hop flows have the same airtime. 
The multi-hop flow in class 1 is allocated a smaller packet size 
(and so less airtime) than the single hop flows. It can also be 



- D1.D2.1 
-D1=1,D2=1e5 
D1.1e5,D2.1 



Number of cells 

Fig. 6. Ratio of airtimes vs. number TV of cells in Parking Lot topology of 
Fig. 5. The y-axis is the ratio of the airtime allocated to the TV-hop flow to 
that allocated to a single hop flow; note that the airtime of the TV-hop flow is 
the sum of allocated airtime in each cell along the flows route. Data is shown 
for three different delay deadline requirements, as indicated in the legend. All 
flows have the same PHY rate. 



class 1 

class 2 

class 3 



w^w, w=54Mbps 

Fig. 7. Coded packet size vs. ratio of PHY rates 103/111 for Parking Lot 
topology of Fig. 5 with TV = 3 cells. Class 1 consists of multi-hop flow 1, 
class 2 consists of single-hop flows 2 and 4, class 3 consists of single-hop 
flow 3; class 1 and 2 flows use PHY rate w bit/sec, the class 3 flow uses a 
PHY rate of W3 bit/sec. All flows have delay deadline D = 1. 



seen that varying the PHY rate for the single-hop flow in class 
3 does not affect the optimal coded packet sizes of flows in 
class 1 and class 2, and hence the airtime of class 1 and class 
2 flows remains the same as W3 is varied. The coded packet 
size of the class 3 flow increases linearly with w 3 /w, and so 
the airtime of the class 4 flow remains invariant as well. 



VIII. Conclusions 

In this paper, we posed a utility fair problem that yields 
the optimum airtime and the coding rate across flows in 
a capacity constrained multi-hop network with delay dead- 
lines. We showed that the problem is highly non-convex. 
Nevertheless, we demonstrate that the global network utility 
optimisation problem can be solved. We obtained the optimum 
airtime/packet size, channel coding rate, and analysed its 
properties. We also analysed some simple networks based on 
the utility optimum framework we proposed. To the best of our 
knowledge, this is the first work on cross-layer optimisation 
that studies optimum coding across flows which are competing 
for network resources and have delay-deadline constraints. 
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Appendix A 
Proof of Lemma 1 



Hence, 



D f-n 



f n f 



tf = p { E E M> 



1=1 

D f n f 



D f n f ~ D f k f 



E E M > D f n f 



1-r/ 



• exp ( -D f n f ■ f ' . ' ) Z 



= exp -D f n f ■ f 



2 

l-r f 



Djnj 

«P \ Vf E 



E 



exp \6fE f [l] 



= exp I —Dfnf ■ Of ^ T * I exp I Dfnf In ( E 



JfE f [l] 



Of^ 1 -In I 2 



= exp -Dfnf 



= exp I -D f n f I Ef[1] [ ^p--6 } 



D,n, 



i= D S n S^P L + 1 



3i 



i-pf 

exp ^—DfUf 



exp ( -D/n/ 
1 + r/ 



1-r/ 



l-f3f 
exp I —Dfnf 



exp ( —Dfnf 
1 + r/ 



ln(l/(l - ft)) 



2/3/ 



In 



1 + r/ 



exp I -D/ri/ 



2(1-/3/) 
'-r/^l-r/ 



1 + r/ Vn^ 1 + r/ 



2 

1-r/ 



• exp ^£>/n/ 

= r ^exp(-D/n/^(B v 2 
exp (-Dfnfl (b (^P-) \\B(p f ) 



e f = 



Appendix B 
Proof of Lemma 2 



E m > 

i=l 
D f n f 



D f n f ~ D f k f 



E E M> D f n f 



2 

1-r/ 



=i 

D f n f 



The binomial coefficients can be bounded as follows: 



1 s$ 



n(n-l)---(n-fc + l) k 
l-2-.-fc 



Appendix C 
Proof of Lemma 3 

For any zeZ, the function f(y, z) is concave in y. Hence, 
for each z, there exists a unique maximum y + (z), which is 
given by 

f(y + (z),z) = max f(y,z) 
=■ 9(z) 

If f(y + {z) 1 z) is a concave function of z, then there exists a 
unique maximiser, which is denoted by z + , i.e., 

z + = argmax f(y + (z),z). 

We show that (y + (z + ), z + ) is an optimum solution to 
Eqn. (9). Since z + is the maximiser of g, we have for any 

z e z, 

g(z + ) > g{z) 

or f(y+(z+),z+) > f(y + (z),z). 

For any given z e Z, y + (z) is the maximiser of f(y, z) over 
all y e y, i.e., 

/(y+(z),z) ^ /(y,z), 
and hence, for all (y, z) e y x Z, 

> f(y + (z),z) > /(j/,z). 

We note that y + (-) maps Z into }>, and hence, (y + (z + ), z + ) G 
y x2. Hence, (y + (z+),z+) is a global maximiser. ■ 
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Appendix D 
Proof of Lemma 4 

Consider the optimisation problem, 
max J2 »/ + ln (! - 2x f(*f)) + ln (! - er ° f cxp(K ^+ 7 ^) 



s.t. 



Ee"/ 
< T c , Vc g C 



We show that the objective function is jointly (strictly) 
concave in (n, I). The objective function is separable in 
(n,f,If), and we show that Xf(If) is convex, and ln(l — 



e --D/exp(n / +/ / )) is concave . 



Since, for Xf £ (/?/, 0.5), If is a monotone function of Xf, 
and If is a monotone function of //, it is clear that // is 
invertible. Note that 



If 


= Hi/) 






dx f 


If 


dx f 


If 


dl f 


9* f (x f ) 


d 2 Xf 


If 


dl) 





0}(x f ) 2 x f (l-x f ) 

Define g(x f ) := ^(1 - x f )6* f (xf) 2 - If. If g(x f ) > 0, 
then Xf(If) is (strictly) convex. Note that g'(xf) = (1 — 
2xf)6*f{xf) 2 + d*f(xf) is increasing with Xf, and hence, 

g(xf) > g(J3f) = o- 

Define h(x,y) = e x+v . consider the function 



Similarly, 



f(x,y) 

df 
dx 
df 
dx 
df 
dy 



\n(l-e- D f 9{x ' v) ) 
D f e- D f g dg 
1 - e- D f3dx 
D } ge- D f g 
1 - e~ D f9 ' 
D f ge- D f 9 
1 - er D s9 ' 



Also, 

d 2 f 
dxdy 



~D 2 g 2 e- 2D fB 



+ 



D f ge 



Dig -D 2 g 2 e- D !3 



(l-e- D f c J) 2 l-e- D f3 
D f9 m ... :•'/' ■/ 



(1 - e- D f3) 2 
D f g(l - e- D f g )e 

~Df9 

(1 - e- D f3) 2 



1 - e -D,g 

[D f ge- 2D f 9 + 

Djg _ (i _ e -D/ff)e-D/fl] 
[D f g-(l-e- D f g )e- D f g ] 



Similarly, one can show that = = dx g y - Define 
i{x) = D f x-(l- e- D f x )e~ D i x . If £{x) > 0, then f(x, y) 
is (strictly) convex. Note that £'(x) = (Df - D f e^ 2D f x ) + 
(D f e- D f x -D f e- 2D t x ) > 0. Therefore, t{x) > £(0) = 0. ■ 

Appendix E 
Proof of Lemma 5 

From Eqn. (21), it is clear that even for a single cell, because 
of the non-zero second term in the LHS, the air-time of flow 
/ given by ^j- is not the same for all the flows /. ■ 
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